Correlating radii and electric monopole transitions of atomic nuclei 
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A systematic analysis of the spherical-to-deformed shape phase transition in even-even rare-earth 
nuclei from 58 Ce to 74 W is carried out in the framework of the interacting boson model. These 
results are then used to calculate nuclear radii and electric monopole (E0) transitions with the same 
effective operator. The influence of the hexadecapole degree of freedom (g boson) on the correlation 
between radii and E0 transitions thus established, is discussed. 
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Electric monopole (E0) transitions between nuclear 
levels proceed mainly by internal conversion with no 
transfer of angular momentum to the ejected electron. 
For transition energies greater than 2m e c 2 1 electron- 
positron pair creation is also possible; two-photon emis- 
sion is possible at all energies but extremely improbable. 
The total probability for a transition between initial and 
final states |i) and |f) can be separated into an electronic 
and a nuclear factor. P = flp 2 , where the nuclear factor 

/>is PI 

P =i 

with R — roA 1 / 3 (tq — 1.2 fm) and where the summation 
runs over the Z protons in the nucleus. The coefficient a 
depends on the assumed nuclear charge distribution but 
in any reasonable case it is smaller than 0.1 and can be 
neglected if the leading term is not too small [l| . 
The charge radius of a state |s) is given by 

<r 2 > S = |<s|X> 2 |s}. (2) 
P =i 

It is found experimentally that the addition of neutrons 
produces a change in the nuclear charge distribution, an 
effect which can be parametrized by means of neutron 
and proton effective charges e n and e p in the charge ra- 
dius operator T{r 2 ). This leads to the following general- 
ization of Eq. (J21): 

1 A 

(r 2 ) s = (s|T(r 2 )|s) = — (s| ^ e fc r 2 |s), (3) 
e n N + e p Z J 

where the sum is over all nucleons and = e n (e p ) if k 
is a neutron (proton). 

An obvious connection between p and the nuclear 
charge radius is established in the approximation a = 
(which henceforth will be made) . Again because of the 



polarization effect of the neutrons, one introduces an E0 
operator of the form [2| 

A 

T(E0) = ]Te fe r 2 . (4) 

k=l 

The p defined in Eq. (|TJ) with a = is then given by p — 
(f |f(E0)|i) /eR 2 . The basic hypothesis of this Letter is to 
assume that the effective nucleon charges in the charge 
radius and E0 transition operators are the same. If this 
is so, comparison of Eqs. and ^) leads to the relation 

f(m) = (e n N + e p Z)f{r 2 ). (5) 

At present, a quantitative test of the correlations be- 
tween radii and E0 transitions implied by ([5]) cannot be 
obtained in the context of the nuclear shell model. The 
main reason is that E0 transitions between states in a sin- 
gle harmonic-oscillator shell vanish identically [l| and a 
non-zero E0 matrix element is obtained only if valence 
nucleons are allowed to occupy at least two oscillator 
shells. This renders the shell-model calculation compu- 
tationally challenging (if not impossible), certainly in the 
heavier nuclei which are considered here. We have there- 
fore chosen to test the implied correlations in the con- 
text of a simpler approach, namely the interacting boson 
model (IBM) of atomic nuclei 0] . In this model low-lying 
collective excitations of nuclei are described in terms of 
iVb bosons distributed over an s and a d (and sometimes 
a g) level which can be thought of as correlated pairs 
of nucleons occupying valence shell-model orbits coupled 
to angular momentum zero and two (and four), respec- 
tively. The number of bosons N\, is thus half the number 
of nucleons in the valence shell. 

As the charge radius operator in the IBM we take the 
following scalar expression in terms of the algebra's [U(6) 
or U(15)] generators 0, Q: 

f{r 2 ) = (r 2 ) c + a N h + ^{r ! h d + 1 n g ), (6) 

where (r 2 ) c is the charge radius of the core nucleus and 
ft*d (■n.p) is the d(g)-boson number operator; a, 77, and 7 
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are parameters with units of length 2 . Then, in analogy 
with Eq. ([5]), the appropriate form of the EO transition 
operator is 

f (EO) = £pJV + *> Z ( V h d + 1 n g ). (7) 

Note that for EO transitions the initial and final states 
are different and neither the constant (r 2 ) c nor iVb con- 
tribute to the transition, so they can be omitted from the 
EO operator. The terms fid and h g in Eq. ([6]) stand for 
the contribution from the quadrupole and hexadecapole 
deformations to the nuclear radius. In a first approxima- 
tion the term in fi g will be omitted from Eqs. (|6l7p . Sub- 
sequently, the influence of the g boson will be explored 
in sdg-lBM, 

Although the coefficients a and r\ are treated as pa- 
rameters and fitted to data on radii, it is important to 
understand their physical relevance. The second term 
in Eq. ([6]) depends linearly on particle number and for a 
not too large range of nuclei it can be associated with the 
isotope shift originating from the average nuclear charge 
distribution which varies as (r 2 ) av « 3r 2 A 2 / 3 /5 @, with 
r = 1.2 fm. An estimate of a follows from 

l«l « |rg ((-A + 2)»/» _ « grgA-Vs, (8 ) 

which for the nuclei considered here (A ~ 150) gives \a\ ~ 
0.2 fm 2 . The contribution of the quadrupole deformation 
to the radius can be estimated as (r 2 ) dcf s» 5/3 2 (r 2 ) av /47r, 
where 8 is the quadrupole deformation parameter of the 
geometric model 6] . An estimate of rj can be obtained by 
associating (r 2 )d f with the expectation value r)(hd)/Nb 
in the ground state of sd-IBM. In a coherent-state ap- 
proximation this leads to the relation 

^TTW w t r ' N ' A ~ 4/ ^ 2 ' ^ 

where use has been made of the approximate correspon- 
dence 8 s» (4Nb/3A)y/nf3 between the quadrupole de- 
formations 8 and 8 in the geometric model and in the 
IBM, respectively [7|. For typical values of N\> ~ 10 and 
A ~ 150 this gives a range of possible rj values between 
0.25 and 0.75 fm 2 . The preceding analysis also reveals 
the dependence of rj on the ratio of the valence to total 
number of nucleons, which is expected due the valence 
character of the IBM. 

To test the correlation implied by (0, we have car- 
ried out a systematic analysis of even-even nuclei in the 
rare-earth region from Z = 58 to Z — 74. Isotope series 
in this region vary from spherical to deformed shapes, 
displaying a more or less sudden shape phase transition. 
Such nuclear behavior can be parametrized in terms of 
the standard sd-IBM Hamiltonian 0] • The details of this 
calculation will be reported in a longer publication Q. 
Suffice it to say here that the procedure is closely related 
to the one followed by Garcfa-Ramos et al. [9[ and yields 
a root-mean square deviation for an entire isotopic chain 
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FIG. 1: Observed (symbols) and calculated (lines) isotope 
shifts (in fm 2 ) for isotopic chains in the rare-earth region. 
The quantity 0.1(74 - Z) fm 2 is added to A(r 2 ) for display 
purposes. Data are from Ref. [rf| . 

which is typically of the order of 100 keV. In deformed 
nuclei one adjusts the IBM Hamiltonian to the observed 
ground-state band and 8- and 7-vibrational bands, if 
known. Care should be taken, however, with the identifi- 
cation of the /3-vibrational band. Our criterion has been 
to select the band with the strongest p 2 (not necessarily 
the first-excited K v = + band) since this is one of the 
characteristic features of a 8- vibrational band [10( ■ 

With wave functions fixed from the energy spectrum 
one can now compute nuclear radii using the operator 
In the study of nuclear phase-transitional behavior it is 
more relevant to plot, instead of the radii (r 2 ) them- 
selves, the isotope shifts A(r 2 }^ = ( r 2 )( A + 2 ) - (r 2 }<- A \ 
The parameter \a\ is adjusted to the different isotopic 
chains (\a\ = 0.22, 0.24, 0.26, 0.13, 0.15, 0.15, 0.11, 0.10, 
and 0.11 fm 2 for even Z between 58 and 74) but other- 
wise all calculated isotope shifts in Fig. [1] are obtained 
with a single parameter rj — 0.5 fm 2 . This value is rather 
reliably determined from the peaks in the isotope shifts 
in the spherical-deformed transitional region and repro- 
duces the essential features of charge radii in all isotopic 
chains with the exception of the light Yb isotopes. 

A more direct way to fix the parameter rj is from isomer 
shifts since this quantity is directly proportional to 11. Of 
the few isomer shifts known in the rare-earth region, only 
the Gd isomer shifts have been measured with different 
techniques to give consistent values and these agree with 
the theoretical results calculated with 77 = 0.5 fm 2 [8|. 

With the parameter rj determined from isotope and 
isomer shifts, one can now compute p 2 values in the rare- 
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TABLE I: Calculated and experimental 0, [T^] p 2 values be- 
tween levels with angular momentum J and initial and final 
energies Ei and E{ (in keV). 



Isotope 


Transition 


J 


p 2 


x 10 3 




£j{ 


Calc 


Expt 


lb,, Sm 


140 


—* 





8 


18 2 




1040 


O Q A 

— > 334 


2 


18 


100 40 


152 Sm 


685 


— > 


o 


57 


51 5 




Oil 


— > Izz 


2 


45 


69 6 




lU2o 


— > ODD 


4 


32 


88 14 




lUoo 


— > u 





3 


0.7 4 




lUoo 


— > OoO 





52 


22 9 


152 Gd 


DIO 


— ► u 





76 


63 14 






— > 344 


2 


86 


35 3 


154 Gd 




— > n 





95 


89 17 




815 


— * 123 


2 


74 


74 9 


156 Gd 


1049 


— > 





58 


42 20 




1129 


— » 89 


2 


46 


55 5 


158 Gd 


1452 


— > 





34 


35 12 




1517 


_> 79 


2 


30 


17 3 


158 Dy 


1086 


— » 99 


2 


47 


27 12 


160 Dy 


1350 


— > 87 


2 


32 


17 4 


i62 Er 


1171 


-> 102 


2 


43 


630 460 


164 Er 


1484 


-> 91 


2 


29 


90 50 


i66 Er 


1460 








15 


2 1 


i70 Yb 


1229 


-> 





36 


27 5 


172y b 


1405 








34 


0.20 3 


174 ffi 


900 


-> 91 


2 


36 


27 13 


176 Hf 


1227 


-> 89 


2 


17 


52 9 


178 Hf 


1496 


-> 93 


2 


36 


14 3 


!82 W 


1257 


100 


2 


51 


3.5 3 


184 w 


1121 


-c 111 


2 


58 


2.6 5 



earth region using the E0 transition operator ([7]). The 
results are compared with the available data in Table [U 

Table Q] illustrates the successes and failures of the 
present approach. In the Sm, Gd, and Dy isotopes, 
the model reproduces the correct order of magnitude of 
the p 2 values for a reasonable choice of effective charges, 
e n = 0.5e and e p = e. In the heavier nuclei, however, 
one observes some glaring discrepancies, most notably 
in 166 Er, 172 Yb, and 182 - 184 w. One possible reason is 
that the strong p 2 to the /3-vibrational band has not 
yet been identified experimentally in these nuclei. This 
seems to be the case in 166 Er where recently Wimmer et 
al. [H[ measured a p 2 value of 127 (60) 1CT 3 to an ex- 
cited + state at 1934 keV. At present this observation is 
difficult to place in the systematics of the /3-vibrational 
band in the Er isotopes and therefore a re-measurement 
of E0 properties seems in order before attempting a re- 
interpretation of these nuclei. Many p 2 values have been 
measured in 172 Yb but none of them is large leading to 
a confusing situation that also deserves to be revisited 
experimentally. Only in the W isotopes it seems certain 
that the observed E0 strength is consistently an order 
of magnitude smaller than the calculation. It is known 
that these nuclei are in a region of hexadecapole defor- 
mation [l4j]- This may offer a qualitative explanation of 



TABLE II: Matrix elements in the classical limit of SU(3) in 
sd-IBM and sdg-IBM. 
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the suppression of the E0 strength, as is argued below. 

From the preceding analysis the following picture 
emerges. Isotopic chains in the rare-earth region exhibit 
a spherical-to-deformed evolution which, at the phase- 
transitional point, is characterized by a peak in the iso- 
tope shifts. Parallel with this behavior there should be 
a peak in the E0 strength from the ground to an excited 
+ state. As emphasized by von Brentano et al. [l5[, an 
inescapable prediction of the sd-IBM is that sizable E0 
strength should be observed in all deformed nuclei. We 
now investigate to what extent this conclusion is 'robust' 
by studying a well-known extension of se?-IBM through 
the introduction of a g boson. 

For a review of studies with the sdg-IBM we refer the 
reader to Devi and Kota [l6j]. A spherical-to-deformed 
transition occurs between the limits U(5) £g> U(9) and 
SU(3) of the sdg-V(15) model [13 ■ Up to a scale fac- 
tor, irrelevant for the subsequent discussion, a schematic 
Hamiltonian, transitional between the two limits, is of 
the form 

H = (l-t)(n d + \h g )--£-Q-Q, (10) 

41V b 

where Q u is the SU(3) quadrupole operator of the sdg- 
IBM [H. The U(5) <8> U(9) limit is obtained for C = 
whereas the SU(3) limit corresponds to f = 1. By varying 
£ from to 1 one will thus cross the critical point £ c w 0.5 
at which the spherical-to-deformed transition occurs. 

Analytic expressions can be derived for the matrix el- 
ements of the operators h s , fid, and n„ for the limiting 
values of £ in the Hamiltonian (fTU|) 8]. It is instruc- 
tive to compare these results to the corresponding ones 
in sd-IBM which is done in Table |TT] in the classical limit 
Ab — > oo of SU(3). Considering first the expectation val- 
ues of h\ in the ground state, we note that d bosons are 
dominant in 0^~ both in sd- and sdg-IBM and that the 
contribution of g bosons in sdg-IBM is small. One there- 
fore does not expect a significant effect of the g boson 
on the nuclear radius, and this should be even more so 
away from the SU(3) limit for a realistic choice of boson 
energies, < ea < e g . Turning to the (0^~ |0^~) matrix 
elements, we see from Table [TT1 that the matrix elements 
of fid and fig in sdg-IBM are of comparable size and differ- 
ent sign. Therefore, while changes in the nuclear radius 
due to the g boson are expected to be small, one cannot 
rule out its significant impact on the p 2 (0t — > 0^ ) values 



4 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 

e e 



FIG. 2: The matrix elements (a) (Of |n ; |0+) and (b) 
{0^|n;|0^) 2 for I — 2 and I = 4 in the spherical-to-deformed 
transition of sd-IBM and sdg-IBM. In sd-IBM the transition 
is from U(5) to SU(3) and sdg-IBM from U(5)®U(9) to SU(3) 
with A = 1.5. The number of bosons is Nb = 8. 

in deformed nuclei. 

This argument can be made more quantitative by 
studying the spherical-to-deformed shape transition of 
the Hamiltonian (jTTJjl . Using a numerical code the 
matrix elements of fid and fi g can be calculated for arbi- 
trary £. For the ratio of boson energies we choose A = 1.5. 
The results are shown in Fig. [5] and compared with the 
matrix elements of fid calculated for the U(5)-to-SU(3) 
transition in sd-IBM. Panel (a) of the figure confirms 
the dominance of the d boson in the ground state of de- 



formed nuclei both in sd- and sdg-IBM. Moreover, the 
expectation value of fid varies with £ in very much the 
same way in both models. In sd-IBM as well as in sdg- 
IBM the sharp increase in (O^ln^O^) 2 is observed around 
Cc ~ 0.5. Up to that point there is essentially no contri- 
bution to p 2 (0^ — > 0i ) from the g boson. Consequently, 
all sd-IBM E0 results up to the phase-transitional point 
are not modified significantly by the g boson. As can be 
seen from Fig.[2L, in the deformed regime this is no longer 
true since, in sdg-IBM, a sharp decrease of (O^ln^O^) 2 
occurs at ( ps 0.5 and (Of |n- 3 |0^) 2 rapidly increases after 
C » 0.6 and dominates (0f\n d \0^) 2 for ( > 0.7. 

In summary, the correlation between radii and E0 tran- 
sitions has been investigated in transitional nuclei. A 
quantitative analysis of radii and p 2 values in the rare- 
earth nuclei seems to validate the explanation of E0 
strength which is based on a geometric picture of the 
nucleus as advocated a long time ago [lfj • The proposed 
correlation depends on the effective charges appearing in 
the radius and E0 operators, assumed to be the same, 
and establishes a method to determine these charges em- 
pirically. The hexadecapole deformation which has only 
a weak effect on the nuclear radius but possibly a strong 
one on E0 transitions, might perturb the correlation in 
deformed nuclei. 
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